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Abstract 

We calculate the order 0(a 2 a z s ) interference effect between the gluon fusion and weak 
boson fusion processes allowed at the one-loop level in Higgs boson plus 2 jet production at 
the LHC. The corresponding one-loop amplitudes, which have not been considered in the lit- 
erature so far, are evaluated analytically using dimensional rcgularisation and the necessary 
master integrals with massive propagators are reported. It is discussed in detail how vari- 
ous mechanisms conspire to make this contribution numerically negligible for experimental 
studies at the LHC. 

1 Introduction 

One of the main tasks for the experimental and theoretical programme in connection with the 
CERN LHC is to investigate the mechanism of electro-weak symmetry breaking. Central to this 
study would be the measurement of the couplings of any observed Higgs scalar to the electro- 
weak bosons. This can be performed either by studying the decays H — > ZZ,WW [1,2] with 
contributions from all production channels, or the production process pp — ► Hjj [3-5] through 
weak boson fusion (WBF) [6], as shown in Fig.[TJa), with contributions from all identifiable decay 
channels. The Higgs plus two jet signature also receives contributions from Higgs boson produc- 
tion through gluon-fusion mediated through a top- loop, as illustrated in Fig. [H(b). However, 
the Higgs plus dijet-sample can be biased towards WBF by suppressing the gluon-fusion chan- 
nel through a combination of cuts, requiring both well-separated jets (effectively suppressing 
the largest component of gluon-initiated processes) and suppressing events with further central 
jets (produced predominantly in the gluon-fusion channel, since here the two tagging jets are 
colour-connected) . 

For the gluon fusion process, the first radiative corrections have been calculated within 
QCD [7, 8] using the heavy top mass effective Lagrangian [9-11]. For the WBF, both the 
radiative corrections within QCD [11-14] and the electro-weak sector [15] have been calculated. 
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Figure 1: (a) The WBF process for Higgs production in the Standard Model and (b) the 
equivalent gluon- fusion diagram mediated through a top-loop. 

The radiative corrections to the WBF channel are small, 3% — 6%, and there is even partial 
numerical cancellation between the QCD and electro-weak contributions. It would therefore 
seem that the Higgs coupling to electro-weak bosons can be very cleanly studied with a Hjj- 
sample. 

Until recently, the irreducible contamination in the extraction of the ZZH-coupling from 
interference between the gluon fusion and WBF processes was ignored in the literature. At 
tree level, such interference is only allowed in amplitudes where the two quarks have the same 
flavours, but their contribution is kinematically suppressed by the requirement of a t <-> u- 
channel crossing, as discussed in Ref. [16]. These interference terms were later also included in 
the calculation reported in Ref. [15], where the full electroweak corrections have been calculated, 
and which also took into account other crossing-suppressed one- loop amplitudes. 

In the present paper we will report on the calculation of the processes allowed at the one- 
loop level which do not suffer from the suppression stemming from the requirement of a t <-> it- 
crossing. As will be explained below, one finds at order C(a 2 af) an interference term between 
the gluon- and Z-induced amplitude which is not allowed at 0(a 2 a 2 s ) by colour conservation. 
The VF-induced amplitudes are crossing-suppressed and therefore not taken into account. The 
diagrams where the vector boson is in the s-channel can be safely neglected because they are 
strongly suppressed by the WBF cuts. 

Given that electroweak corrections to the WBF process, which are formally an order 0{a) 
correction to an 0(a 4 ) process, have been shown to be relevant for this important process [15], 
simple power counting alone suggests that the size of the irreducible contamination due to the 
discussed interference effect should be checked. We will elaborate below that arguments in the 
literature which are based on simplified assumptions do not capture all effects found by doing 
the full one-loop calculation. 

In the following section we will briefly sketch the calculation before discussing our results in 
section [3l which are summarized in the conclusions. The appendix contains an extensive list of 
the master integrals needed for this calculation. Most of these integrals have not been reported 
in the literature so far. 
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Figure 2: Example of contributing one-loop interference terms: (a) M. g z-M* g and (b) M. gg Ai* z . 
There are four contributing topologies for each gluon-fusion and Z-fusion process. 




Figure 3: One of the real emission processes which contributes at the matrix element squared 
level. 

2 The Calculation 

We set out to calculate the helicity amplitudes necessary to form the loop interference terms 
and the real emission contributions. Sample diagrams are shown in Figs. [2] and [3l As discussed 
in Ref. [16] this is the lowest order contribution to the interference between ZZH and ggH- 
processes for non-identical quark flavours and helicity configurations, and for identical quark 
and helicity configurations the loop amplitudes are the first order which does not require a 
kinematically disfavoured crossing. 

The amplitudes have four non-zero helicity components, which we label by + + +- 1-, , 

— h H — and H h Due to parity invariance of the kinematical part of the amplitudes, only 

two of them are independent. By using the spinor helicity formalism we have defined projection 
operators on each of these amplitudes. In practice, we calculate the amplitudes with all momenta 
incoming (and therefore summing to zero); to map to physical scattering kinematics, crossing 
relations are applied easily in the end. 

The leading order amplitudes, denoted by Mz and M. g (see Fig. [2]), are proportional to a 
colour singlet and a colour octet term. The colour singlet is formally of order 0(a 2 ) whereas 
the octet is of order 0{a 2 s ). The virtual corrections, which we call M. g z and M. gg respectively, 



are mixtures of octet and singlet terms. For the intereference term we need to consider only the 
octet part of Mgz and the singlet part of M gg - One finds that only four one-loop five-point 
topologies for each amplitude survive this colour projection. As was already pointed out in 
Ref. [16], the colour singlet cannot interfere with the colour octet tree amplitude for different 
quark flavours. However, a new colour channel opens up at order 0(a 2 a^) which is neither 
flavour nor kinematically suppressed. 

The loop amplitudes require the evaluation of one-loop five-point tensor integrals with partly 
massive propagators and external legs. We apply the reduction algorithm outlined in Ref. [17,18] 
to express each Feynman diagram as a linear combination of 1-, 2-, and 3-point functions in 
D = 4 — 2e dimensions and 4-point functions in D=6. The same algorithm has been successfully 
applied to a number of one-loop computations and further details can be found elsewhere [19-22]. 
The coefficient of each integral, which is a rational polynomial in terms of Mandelstam variables 
s ij = {Pi +Pj) 2 and masses, was evaluated symbolically using FORM [23] and simplified using 
Mathematica [24]. Both steps were fully automated. The algebraic expressions were checked 
by independent implementations, both amongst the authors and with another group [25]. 

After the algebraic reduction, all helicity amplitudes for both cases, gluon and weak boson 
fusion, were obtained as linear combinations of a certain number of scalar integrals. We choose 
this basis of so-called master integrals (Mis) in accordance with Ref. [26], i.e. our Mi's are 
D-dimensional two-point and three-point functions (I®,!®), and (Z)+2)-dimensional four-point 
functions (I^ +2 ). Schematically 

M = Y, W,(K,ma}) , Ij € {I?, if, if +2 } (1) 

where the summation over a indicates the summation over different argument lists {s a ,m a } of 
the relevant MI. The conventions for the arguments and the analytic forms of these are given 
in the appendix. No one-point functions appear in the reduction, and also two-point functions 
are absent in the amplitudes of M g z- Furthermore, coefficients of some of the integrals which 
arise in several topologies sum to zero: if the tree resulting from a cut of an internal line of a 
master integral corresponds to helicity forbidden tree level processes one can immediately infer 
the vanishing of the corresponding coefficient. In our algebraic tensor reduction approach we are 
able to verify such cancellations and enforce them analytically before the numerical evaluation 
of the cross section (for a nontrivial example see Ref. [26]). 

The coefficients we obtain through this procedure can be too large to be of use in a printed 
form (their simple polynomial structure means that numerical evaluation is, however, fast). We 
performed several algebraic checks of relations between coefficients of different topologies and 
helicity configurations. The coefficients are included as supplementary material to this paper. 

As most of the required integrals are not provided in the literature, we have evaluated 
representations in terms of analytic functions valid in all kinematic regions. We give our result 
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> (600 GeV) 2 





Table 1: The cuts used in the following analysis which bias the Higgs Boson plus dijet sample 
towards WBF. The suffices a, b label the tagged jets. 

for these integrals in the appendix, as they might be of use for other calculation^. 

The IR structure of the interference term is very simple. It is easily extracted from the result 
by focusing on the IR divergent triangles. All Mis with single poles drop out. Only triangles 
with double poles survive. This results in an expression 

~ a s C F [(-s 13 )^ + (-s 24 )- e - (si2)- e ~ (-s 3 4)- £ ]/e 2 

in which all double poles cancel and only sub-leading soft divergences survive. 

The virtual corrections to the interference term have to be combined with the real emission 
part shown in Fig. In accordance with the virtual corrections, the collinear IR divergences 
from the three-parton final states integrate to zero, leaving only a soft divergence proportional 
to 1/e in dimensional regularisation. Due to the simple structure of the divergences, we have 
used the phase space slicing method [28,29] to isolate the IR divergences from the real radiation 
part. We have checked that the remaining single poles cancel exactly when combining the real 
emission with the virtual part. The phase space integration and the numerical evaluation of 
integrals and coefficients is coded in a C++ program allowing for a flexible implementation of 
cuts and observables. 



3 Results 

This study aims at investigating a possible pollution of the clean extraction of the ZZH vertex 
structure by the interference terms. Therefore we will apply the cuts summarised in Table [TJ 
which are generally used for the selection of WBF events [30] over the gluon fusion "background" . 
Our input parameters for the numerical studies are taken from either the parton density function- 
fit [31] in the case of as(M^) and the Review of Particle Physics [32] for the others. 

a s {M 2 z ) = 0.1205 , g 2 = ° F 1 2 , G F = 1.16637 x 10~ 5 GeV" 2 

V2 8M W (2) 

M z = 91.1876 GeV , M w = 80.425 GeV , sin 2 9 W = 0.2312 . 



1 Ol course all finite integrals can in principle be evaluated by using the LoopTools package [27]. The D = 6 
boxes can be written as linear combinations of 3- and 4-point functions in D=4. If IR divergences are present, a 
small regulator mass must be used, but its dependence can be made arbitrarily small as the D=6 box integral is 
IR finite. 
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Figure 4: The A^-distribution for different helicity-configurations of the valence quarks only. 
The purple histogram labelled "Sum" indicates the sum over the four contributions shown. The 
sum over all valence quark flavour and helicity assignments is shown in the black histogram. 

We have checked that variations of the numerical values chosen for the WBF cuts have no impact 
on the conclusions, nor does the exact value of the Higgs boson mass, which we set to 115 GeV 
unless otherwise stated. The same is true for the choice of parton sets; we choose to use the 
NLO set from Ref. [31], and use 2-loop running for a s , in accordance with the chosen pdfs. 

We observe that in all the flavour and helicity channels, the finite contribution from the 
3-parton final state is numerically negligible. In fact, in the case at hand, the only role of this 
real emission is to cancel the divergences which arise from the one- loop diagrams. 

As an interference effect proportional to 2Re(MggM* z +M g z-Mg), the result is not necessar- 
ily positive definite. In fact, the sign of the interference contribution depends on the azimuthal 
angle between the two tagging jets, A<f>jj. Because of the event topology with two well sep- 
arated jets, it becomes possible to define an orientation of the azimuthal angle which allows 
observability in the whole range of [— 7r, 7t), as pioneered in Ref. [6,33]. ^(pjj is then defined 
through 

\P+T I \P~T I COS A( fe = P+T ■ P-T > ( 3 ) 

2\p +T \\p- T \smA(j)jj = e^ vpcT b^p L lb p _p'L 1 

where 6+ (&_) are unit vectors in positive (negative) beam direction, and likewise for the jet 
momenta p±. The cuts ensure that the two tagging jets lie in opposite hemispheres. 

Figure [5] displays the contribution to the distribution in Acfijj from the interference terms for 
various helicity and flavour configurations of the valence quarks only, for a Higgs boson mass of 
115 GeV. Figure [5] displays the contribution to the distribution in A(f>jj including sea quarks. 

Due to the oscillatory behaviour, the total integrated cross section does not at all tell the 
full story about the size of the impact on the A (^-distribution; for example, the integral of 
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Figure 5: The A<^, -distribution for various flavour and helicity-configurations. The purple 
histogram labelled "Sum" indicates the sum over the four contributions shown. The sum over 
all flavour and helicity assignments including all sea flavours is shown in the black histogram. 

the contribution from the sea and valence up-quarks with negative helicity, denoted by u~u~, 
is +5 ab, while the distribution peaks at more than 150 ab/rad. If only valence quarks are 
considered, the integral is —30 ab, while the distribution peaks at ~ 90 ab/rad. 

There is an accidental cancellation of sea and valence quark contributions which leads to the 
fact that the sum over all flavour and helicity assignments peaks at around 2 ab/rad only, with 
an integrated effect of 1.19 ± 0.07 ab, where the error is due to the numerical integration. Note 
that the integral of the absolute value of the (f>jj distribution, 




is a useful measure of the impact of the interference effect on the extraction of the ZZH-vertex. 
This integral evaluates to 9.1 ± 0.1 ab, an order of magnitude larger. The total integral over the 
the absolute value of the fully differential cross section leads to 29.59 ± 0.07 ab. 

As can be readily seen, there is a cancellation between the contribution from each flavour and 
helicity assignment; this is because the sign of quark couplings to the Z-boson becomes relevant, 
since it is not squared for the interference. The flavour- and helicity sum for each quark line 
therefore leads to some cancellation, resulting from the weak charges listed in Table [2j This 
cancellation was discussed in Ref. [34] for the valence content where it was pointed out that if 
sin 2 6w = 1/4 and the it- valence content of the proton was exactly twice the d- valence content 
for all Bjorken x, the contributions from uu, ud and dd would sum to zero. However, neither of 
these approximations is exactly true: sin 2 9\y = 0.2312 and the parton distribution functions are 
shown in Fig. [6l The cancellation stemming from the SU{2) x [/(l)-flavour and helicity sums 
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-2/3 
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-1/2 


1/3 
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1/3 



Table 2: SU{2) x £/(l) charges for the Zgq-couplings, [a/^ + &/,a sin 2 (# w )] 

of the valence quarks can be effectively studied by calculating the ratio 

Z^2 c f,\( x iQ 2 )//~2 \ c fA x 'Q 2 )\ where 

(4) 

Cf,x(x,Q 2 ) = f{x,Q 2 ) [a fi x + bf tX sm 2 (9 w )], 

where f(x, Q 2 ) is the relevant parton distribution function. This ratio is plotted in Fig. [7] (left). 
This figure also shows the ratio of the u-valence to the <i-valence pdf, which influences the 
cancellation. It follows that in QCD the cancellation due to SU(2) x U (l)-flavour and helicity 
traces amounts to roughly 10 _1 in the most relevant regions of the pdfsj. If one includes the 
sea quarks in this argument, one sees that the ratio Eq. (j4]) is modified significantly, see Fig. [7] 
(right). Its change of sign in the relavant x-region leads to a further reduction of the interference 
term after integration. 

Using the same cuts and value for the mass of the Higgs boson as in the present study, we 
have checked that the total contribution to the A</>,-,-distribution from the leading order WBF 
process (both Z and W + ^~ included) is relatively flat at around 240 fb/rad. Therefore, the 
result of the interference effect reported here is unlikely to be measurable. 

The smallness of the overall effect is in fact also a result of the complex phases arising from 
the full one-loop calculation of the amplitudes. To illustrate this, we calculate the average value 
of 

\Re{M gZ M g + M gg M z )\ 
\M gZ M* g \ + \M gg M* z \ • 1 ' 

The average over phase space for this quantity is roughly 20%, which illustrates that the relevant 
products and sums for the interference effect project out only a small component of the full 
complex loop amplitudes. 

We have checked that none of the sources of suppression discussed above depends severely 
on the mass of the Higgs boson; since the amplitudes themselves depend on this parameter 
only weakly, the effect of increasing the Higgs boson mass is basically nothing but a reduction 

2 We were not able to reproduce the one order of magnitude larger suppression factor reported in Ref. [34] 
within the naive Quark Model. 
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Figure 6: The required parton distribution functions [31] in the relevant region of Bjorken x for 
Q 2 = 400 GeV 2. 

of available phase space. This reduction however is very small, since the partonic centre of 
mass energy is dominated by the contribution from the jets rather than the Higgs boson. The 
exact numerical value chosen for the cuts also does not affect the relative importance of the 
interference effect — the effect on the interference and the WBF signal is similar, so the relative 
importance of the interference is largely unchanged. 

For completeness we list in Table [3] the integral of the absolute A (^-distribution for various 
choices of the renormalisation and factorisation scales. 

We chose the factorisation and renormalisation scales as in accordance with the natural scales 
in the relevant high energy limit (as in Ref. [35]), i.e. the factorisation scales are set equal to the 







Integral of \da/dA(/)jj [ab] 


Pa T 


Pb T 


9.1 ±0.1 


m H /2 


m H /2 


13.9 ±0.6 


mn 




9.2 ±0.4 


2m,H 


2l7lH 


6.3 ±0.3 



Table 3: The dependence of the interference effect on the choices for factorisation and renormal- 
isation scales, with a fixed mass of the Higgs boson of 115 GeV. The variation is slightly larger 
than what would be expected from renormalisation scale variations only, since the variations in 
factorisation scale impacts the cancellations between SU(2) x [/(l)-charges by slightly altering 
the ratio defined above Eq. ([!]). 
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Figure 7: The pdf-weighted sum of the SU(2) x L r (l)-charges divided by the sum of the absolute 
weighted charges for valence quarks only (left) and valence and sea quarks (right). The left figure 
illustrates how the apparent almost complete cancellation expected from the Quark Model is in 
fact not as severe once the pdfs are taken into account, as discussed in the text. The insert shows 
the ratio of the valence u-quark distribution to the valence d-quark distribution for the pdf set 
MRST 2004 (NLO) [31]. The right plot shows that the inclusion of sea quarks does actually 
alter the pdf-weighted sum in the relevant x range such that there are additional compensations 
when integrating over x. 

transverse momenta of the relevant jet, and the renormalisation scale for the strong couplings 
are chosen correspondingly, i.e. one a s evaluated at each value of the transverse momentum of 
the jets, and one at the Higgs mass. However, varying these has no impact on the conclusions. 

4 Conclusions 

We have presented the calculation of the loop-induced 0{a 2 a^) interference effect between the 
gluon fusion and weak boson fusion processes in Higgs boson plus two jet production at the 
LHC. 

In the context of the weak boson fusion cuts we have evaluated all relevant one-loop diagrams 
algebraically and have obtained an analytic representation of the interference term as a linear 
combination of scalar one-loop integrals. The analytic result for all the necessary integrals is 
presented for general kinematics such that it can be used in other computations. 

Our expressions have been coded into a flexible computer program to test speculations in the 
literature about the size of this interference contribution. We do confirm by explicit calculation 
that this contribution is too small to contaminate the extraction of the ZZff-coupling from WBF 
processes. Interestingly the effect which survives comes dominantly from the virtual corrections. 
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We have analysed in detail why this contribution is so small, and instead of a single effect we 
rather find a conspiracy of several mechanisms which can only be completely assessed having 
the full NLO calculation at hand. 
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The mechanisms basically are 

• accidental cancellations between the sea quark and valence quark contributions 

• compensations between different weak isospin flavours of the valence quark contributions 
due to their SU(2) x U{\) couplings in combination with their weights from the (valence) 
quark content of the proton 

• cancellations due to destructive interference of the phases from the different contributions. 

The exact impact of these partly accidental effects has until now not been quantified thoroughly 
and was very hard to assess without an explicit calculation. 

As a final comment we would like to point out that anomalous couplings which affect the 
phases could change the interference pattern substantially. However, the first two cancellation 
mechanisms still being present, the overall contribution is still expected to be undetectably 
small. 
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A Analytic Results for Master Integrals 



The appendix contains the Master Integrals (Mis) which occur in the reduction of the one-loop 
pentagon diagrams encountered in the given calculation. The analytic results of the integrals 
containing only massless propagators have appeared in the literature, see e.g. [17,36-38]. Those 
with massive propagators have been calculated for this project. All finite integrals can also be 
calculated by the LoopTools package [27], based on [39-41] which we used for checking purposes. 

The conventions we use for the scalar triangles and boxes listed below are different from the 
ones defined in [17,18,36,38], as we follow the LoopTools conventions for the argument lists, in 
order to comply with the "Les Houches Accord on Master Integrals" . 

To be specific, for a general Appoint integral as shown in Fig. [8j we use 

r d D k n A - D 

tD( r _ i . 2i \ / a K H 1 



I N ({s jl ... jn };{m i }) J i7r o/2 [(k + n)*-r4+iS\...[(k + r N ) 2 -m 2 N + i8\ ^ 

n /• / N \ An v^A r 




2 'J \ LL I I 1 c -K\ N ~T 
[-^ZiOijZj - 10) 

where = Yl)=iPj an ^ $ij = { r i — r j) 2 — m i — m2 - The results for the integrals are charac- 
terised by the invariants Sj 1 ,,j n = {p.j 1 + ... + Pj n ) 2 and m 2 , where the list {sj 1 ...j n } contains 
the invariants defined by ra-particle cuts of the diagram. To map the labelling of 5^ resp. 




Figure 8: Momentum and mass assignments for a general Appoint one-loop graph. 

Fig. [8J to the LoopTools conventions, the argument lists of triangles and boxes are given by 
Is(s 1 ,S2, s 3 ; raj, m\ , m|) and if (si, s 2 , s 3 , s 4 ; s 12 , s 2 3;ml, mf, m\,m 2 z ). 

The integrals listed below will also be posted shortly to the LoopForge integral database at 



http://www.ippp.dur.ac.uk/LoopForge, together with some numerical benchmark points. The 
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triangles (|A3|) . (|X6|) . (lATHl) . ([Aljl) and the boxes (lAlil) . (|Al5l) . fAlB]) . (lAl8|) . (lAl9j) can also 
be found at http://qcdloop.fnal.gov/ for certain kinematic regions. 

For the analytical representations given below we use the following auxiliary functions. 
The Kallen function: 

X(x, y, z) = x 2 + y 2 + z 2 - 2xy - 2yz - 2zx, (A.2) 

the 77-function: 

r](x,y) = log(xy) - log(x) - log(y), (A.3) 

and the R- function [39,40]: 

1 

log(y - z T i$) ~ log(yo - z T i$) 



R(y ,z±i5) = J dy 



y-yo 



Li 2 (>i) - Li 2 (z 2 ) + m log(zi) - 7] 2 log(z 2 ) (A.4) 



, yo yo-i 

where zi = — , z 2 



Vi = v(-zTi5A/(yo - zTiS)) , ??2 = 7/(1 — z^iS, 1/ (3/0 — 2 T *^))- 

We work in general in = 4 — 2e dimensions but give all formulas only up to O(e ). We 
abbreviate infinitesimal displacements in the Mandelstam variables as s = s + i5. 



Triangle integrals 

In the figures below, a single (double) internal line represents a massless (massive) propagator, 
while a single (double) external leg represents one for which p 2 is zero (non-zero) . 



Triangle if (si, 0,0; 0,0, 0) 




Pi 



Jf(si, 0,0; 0,0,0) 



r(l + e) 



si 



1 1 ^ 

~2 - - fog 



7T 

"6 



21 



(A.5) 
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Triangle I°( Sl , s 2 , 0; 0, 0, 0) 




/ 3 D (si,s 2 ,0; 0,0,0) 
r(l + e) (1 



(A.6) 



si - s 2 I e 



log 



-^2 

M 2 



log 



si 



log^ 



-52 

M 2 



log' 



-si 



Triangle /f(si,s 2 ,3 3 ; 0,0,0) 




/f (si, 5 2 , s 3 ; 0,0,0) 



v / A(si,s 2 ,S3) - ^ «i 



2Li 5 



+ 2Li 2 

y+J 



x + y 3 



4* fen*" 



2 \y+ 



1 



log 2 



(A.7) 



with 



x± 



gl + ^3 - 52 T \/ A(si, S 2 , ^3) - ^ $1 

2 si 



y± = 1 - ^t- 



Note that the permutation symmetry of the integral in si,s 2 ,S3 is preserved, although not 
manifest due to the choice of the denominator in x±. 
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Triangle /f (si, 0, 0; 0, M 2 , 0) 




/ 3 w (si,0,0;0,Af 2 ,0) : 
+Li 2 • h 



si — M 2 / 2 



r(l + e) f 1 



-s\ i e 



log 



-si 



+ M 2 



-gi+M 2 \ 



log 2 



M 2 

7^ 



Triangle lf(si, s 2 , 0; 0, M 2 , 0) 




/^( Sl , S2 ,0;0,M 2 ,0) 



+Li 2 



r(i + e ) ri 



Sl 



•51 



M 2 



s 2 — si (^e 
- Li 2 



•52 



•52 



M 2 



-f-lOE 



-s 2 + M 2 \ 1 



A* 



-si + AP 



log 



Triangle I 3 D (si, s 2 , 0; 0, 0, M 2 ) 




/ 3 B (si,s 2) 0;0,0,M 2 ) 
with rro 



[r{xq, xi)-^-+ Li 2 ( 1 

s 2 - si I 6 



•51 



si - s 2 



, Xl 



Sl 



si - s 2 + M 2 
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Triangle /f (si, s 2 , 0; 0, M 2 , M 2 ) 




/^(si, S2 ,0;0,M 2 ,M 2 ^ 

^{ Li2 (i&)" Li2 (i 



n / m 2 

Li 2 ( — J + i? (x ,— 



+R(1 - x Q , x-) - R(x Q , x-) - r)o log 



1-xq 
-x 



(A.ll) 



with 



8l 1 / / 4M 2 

= 1 , a;± = - 1 ± W 1 — 

S2 2 1 V s 2 



V ( 1 - TTo^O, 



M 2 



M 2 M 2 — s 2 x (l — x o) 



Triangle if (s 1} s 2 , s 3 ; 0, M 2 , M 2 ) 




/f(si,s 2 ,s 3 ;0,M 2 ,M^) 



(A.12) 



i?(x_,xi) — R(x + , x\) + i?(l — x_, xo) — — x+, xo) 



a/A(si,s 2 ,s 3 ) 
-R(x-, x ) + i?(x+, x ) -r?_ log 



1 — x_ 



?7 + log 



1 — Xj 
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with 



x± 



Sl + S 2 ~ S3 =F \ / X(s 1 ,S 2 ,S3) 

si -M 2 



x = 77 1 - A /-I — , Xi 




Sl - S3 



T}± = 7] ( -s 3 x± - Sl(l - x±) + M 2 - Z(5, 



M 2 - x±(l - x±)s 2 



In this formula we assume momentum conservation and at least one positive invariant Sj. The 
latter condition is always guaranteed if the triangle graph is a subgraph of 1 — > n or 2 — ► n 
scattering kinematics. These lead to a positive Kaellen function: A(si,S2,S3) > 0. 



Box integrals 

The results are given here for D=A-2e dimensional boxes. The conversion relation to 6-dimensional 
boxes is achieved by the formula 

4 

I? = E^ T Si + (D~3)B if + 2 , (A.13) 

8=1 

where 6, = X^jWL^iy 1 and -B = ^2%=i -^fi denotes the "pinch" triangle stemming from a box 
where the zth propagator is omitted. All pinch triangle integrals needed by eq. (IA.13j) are given 
above. We use S12 = s and S23 = t in the following. 



Box if (si, 0,0,0; a, t; 0,0,0, 0) 




if (si, 0,0,0; s,i; 0,0, 0,0) 



r(l + e) f 2 2 



si 



e 2 e 



1 7?) + log G? 



log 



-si 



+ 2Li 2 (l-|-j + 21og f^j log f^- 



(A.14) 



+ log 2 



(^) 



A 2tt' 
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Box Jf (si,s 2 , 0,0; s,t; 0,0, 0,0) 



If (si, s 2 , 0,0; 0,0, 0,0) 



r(l+e) ( 1 + 1 



si 



-2Li 2 1 



e 2 e 



si 



log 



l ?)- 21 ° g (? 



_|j +log m lo J| 



' 7? ' ~ ^ lc 



M 2 



Box If (si,0, s 3 , 0;s,t; 0,0, 0,0) 



if (si,0, s 3 , 0;s,i; 0,0, 0,0) 



r(l + e) (2 



st — S1S3 I e 



log 
- log 



-si 



+ log 



-S3 

M 2 



log 



log 



/' 



-2Li 2 I 
+2Li 2 ( 



si 
s 

Si S3 



log 2 



S3 



-S3 



2 Li, 1-^ 



I- j 



st 



log 2 



J +2 , f,f,lo g( l 



2Li 2 (^1 

_ hli) 

st J 
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Box If (si, 0,0,0; s,t;0,M 2 , 0,0) 



We find for the 6-dimensional integral 



where 



lf +2 (si,0,0,0;s,t;0,Af 2 ,0,0) 



-t + M 2 
t(si — s — t) 
M 2 
t(si - s) 



R(x ,x 2 ) + Li 2 ( 1 - — 
i? (xi,x 2 ) + Li2 (l 



351 



x 



S + t — Si 

s 



x = x + 



id 



t-M 2 ' 



Xl 



s — Si 

The integral in D = 4 — 2e is obtained by 



xi = xi — i5 , X2 



7T 

IT 



7T 
~~6 



Sl + M 2 



(A.17) 



if (si, 0,0,0; s,t;0,M 2 , 0,0) = 6i I 3 (s, 0, 0; 0, 0, 0) + b 2 I 3 (si, t, 0; 0, M 2 , 0) (A.18) 
+6 3 / 3 (si, 0, s; 0, M 2 , 0) + 64 h(t, 0, 0; 0, M 2 , 0) + B if +2 (si, 0, 0, 0; s, t; 0, M 2 , 0, 0) 



with 



h 

in 



1 



M 2 - t 
t 



si-t t(s-si) + M 2 (s + 2t-si) 



s (M 2 - t) 
_ 2t(s + t-s 1 ) 
s (M 2 - t) ' " ~ s(M 2 -t) 2 ' 



s (M 2 - t)< 



B 
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Box if (si,s 2 , 0,0; s,t;0,M 2 ,M 2 , 0) 



with 



-^polc 
Ifinitc 



X± 



Xl 



if (si,s 2 , 0,0; s,t;0,M 2 ,M 2 , 0) 



-^polc ~t~ -^finite 

si - (s + i - si)M 2 



r(i + e) 



log( j±^!y l0 j-^+n 



M 2 
■log 2 



2l '-t + M 2 \ , 2 (M 2 \ 2 /-s l + M 2 s 

log 



2 ~ lo * 

a** y V A» / V M' 
2 i?(s 0> 5i) + 2 i?(x , x 2 ) 

R(l — xi, X-) — R(xi,x-) — R(l — X2, X-) + R(x2, X-) 

1 — X\ 



-s + M 2 J 
2 f-h 



+ log" 



+ 



log (-h/fJ 2 ) + log(x+ - xi) + log(xi - x_) 
log (-s 2 //i 2 ) + log(x + - x 2 ) + log(>2 - »-) 



log 
log 



-Xl 

1 - x 2 

-X2 



2 \ V s 2 



S + t — 3\ 



M 2 M 2 -i5 s-M 2 s-M 2 
—— , Xi = , x 2 = , x 2 



t 



s — Si 



s — Si 



Box if (si,0, s 3 , 0; s,t; 0,0, M 2 ,M 2 ): 



2 A/r2\. 



If (si, 0,53,05 8, i; 0,0, AT*, AT 



2 j\ / t2^ 



I- — li 
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with the auxiliary integrals^ 




M 2 



R(x±,x ) 



M 2 



t 



s 



) 



(A.21) 



where 



det(S) 



J = 



[st - 
2(s 3 



sis 3 + M 2 (s - t)} 2 + AM 2 (st - sis 3 )(-s + si + M 2 ) 

- t)[st - Sl s 3 + M 2 (s- t)] + 4(st - sis 3 )(-s + si + M 2 ), 



x± = 



(si 




Checks on the integrals 

We tested the Master Integrals by comparing our results numerically with LoopTools [27]. For 
the IR finite box and triangle integrals this is straightforward. The IR divergent box integrals 
have been checked indirectly by mapping to the 6D case which is IR finite. Using eq. (|A.13|) 
the 1/e poles cancel when combining the D-dimensional box with the triangle pinch integrals. 
The same expression can be evaluated with LoopTools by using a mass regulators in the IR 
divergent integrals. In the given IR finite combinations the cut-off dependence is polynomial and 
can be made arbitrarily small. We have tested the formulae in all kinematically distinguishable 
regions. For example, I^(si, s 2 , s 3 ; 0, M 2 , M 2 ) has been checked in the regions resulting from 
combining the conditions {s\ < 0, < s 1 < M 2 , M 2 < si), (s 2 < 0, < s 2 < 4M 2 , AM 2 < s 2 ) 
and (s 3 < 0, < s 3 < M 2 , M 2 < s 3 ) such that the Kaellen function is positive in line with the 
comment below equation (|A.12p . 
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